arXiv:hep-th/0112030vl 4 Dec 2001 


FINITE TEMPERATURE INDUCED FERMION NUMBER 


Gerald V. Dunne 

Department of Physics, University of Connecticut, Stores, CT 06269-30f6, USA 


The induced fractional fermion number at zero temperature is topological (in the sense that it 
is only sensitive to global asymptotic properties of the background field), and is a sharp observable 
(in the sense that it has vanishing rms fluctuations). In contrast, at finite temperature, it is shown 
to be generically nontopological, and not a sharp observable. 


I. INTRODUCTION 


The phenomenon of induced fermion number arises due to the interaction of fermions with nontrivial topological 
backgrounds (e.g., solitons, vortices, monopoles, skyrmions), and has many applications ranging from polymer physics 
to particle physics @ 1 - At zero temperature, the induced fermion number is a topological quantity, and is related 
to the spectral asymmetry of the relevant Dirac operator, which counts the difference between the number of positive 
and negative energy states in the fermion spectrum. Mathematical results, such as index theorems and Levinson’s 
theorem, imply that the zero temperature induced fermion number is determined by the asymptotic topological 
properties of the background fields |i] [ 9 ). This topological character of the induced fermion number is a key feature 
of its application in certain model field theories. At finite temperature, the situation is very different 00- In this 
talk I explain why at finite T the induced fermion number is generically nontopological, and is not a sharp observable. 

This work has also been partly motivated by recent results that certain anomalies in T = 0 field theory become 
much more subtle at nonzero T. For example ||12|| , in even dim. spacetime, finite T anomalous 7r° decay amplitudes 
are T dependent even though the chiral anomaly (a topological object that is related to the anomalous 7r° decay 
amplitude at T = 0) is independent of T. And, in odd dim. spacetime, the topological Chern-Simons term is the only 
parity-violating term induced in the T = 0 effective action, but at finite T there are infinitely many parity-violating 
terms 13-151. 


II. INDUCED FERMION NUMBER 

The induced fermion number is an expectation value of the fermion number operator, TV = J dx ^[\U, \H]. For a 
given classical static background field configuration, the fermion field operator T can be expanded in a complete set 

of eigenstates of the Dirac Hamiltonian H. The T = 0 fermion number, (TV) 0 = (0|iV|0), is related to the spectral 

asymmetry of the Dirac Hamiltonian: 

1 1 

(TV) 0 = —- / dE a(E) sign(£’) = —-(spectral asymmetry). (1) 

• J — OO " 

Here cr(E) is the spectral function of H\ 

a(E) = - 1m Tr ( * . ) . (2) 

7 r V E — E — le J 


So, (TV )0 essentially counts the number of positive energy states minus the number of negative energy states. 
At nonzero temperature T, the induced fermion number is a thermal expectation value 


(TV) t 


Tr (e-P H TV) 
Tr (e~P H ) 


dE a(E) tanh 



( 3 ) 


where /? = A. This finite T expression (j|) reduces smoothly to the zero T expression (jl;) as (3 —> 00 ; in fact, finite 
temperature provides a physically meaningful and mathematically elegant regularization of the spectral asymmetry. 
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Note that the spectrum of the fermions in the static background has nothing to do with the temperature. All 
information about the fermion spectrum is encoded in the spectral function. Thus, if one knows a(E ), one immediately 
has integral representations for both {N)q and ( N)t■ In fact, one only needs to know the odd part of <r(E), as is 
clear from ([if]) and (jpj). It is convenient to use the explicit form (|j) of the spectral function to write (||) as a contour 
integral involving the resolvent of the Dirac Hamiltonian 


(N) t 


1 / A iv 

2 J c 2m 




(4) 


where C is the contour (—oo + it , +oo + it) and (+oo — it, —oo — it) in the complex energy plane. Thus the resolvent 
is the key to the computation. 


III. TWO 1+1 DIMENSIONAL EXAMPLES 

Before going to the general case, it is instructive to consider two important cases in 1+1 dim field theory. These 
two examples are very similar at T=0, but turn out to be very different for T > 0. Consider an abelian model of 
fermions interacting via scalar and pseudoscalar couplings to two (classical, static) bosonic fields 0 i and cj> 2 : 

£ = i-0^0 — -0 (0! +*75 02)^ (5) 

We distinguish between two different cases: 

(a) kink case jlj: 0 i = m is constant, but 02 (a;) has a kink-like shape (for example 02 ( 2 ;) = tanh(x)), with 
02 (± 00 ) = ± 02 - 

(b) sigma-model case ||: 0 i (x) and 02 (x) are constrained to the “chiral circle”: 0 f + 0 | = to 2 . 

The kink case (a) is relevant to the polymer physics applications, while the sigma-model case (b) is a low-dimensional 
analogue of sigma models used in 3+1 dim. particle physics models [Q. 

In each case, we define the angular field 


9(x) = arctan [02 (a;)/0i (a;)] (6) 

In the sigma model case (b), the angular field 9(x) has the interpretation of a local chiral angle, and is chosen to 
have a kink-like shape, with asymptotic values 0(±oo) = ±6. Note that in the kink case (a), the angular field 
9{x) = arctan[02 {x)/m\ also has a kink-like shape. 

For both the kink case and the sigma-model case, at zero T, the induced fermion number is 

( N )o = -^f dxd '( x ) = ^ ( 7 ) 

which is topological in the sense that it only depends on the asymptotic values of 9{x), not on its detailed shape. In 
the kink case, this calculation is easy because the even part of the corresponding Dirac resolvent is known exactly 

=-( 8 ) 

e (to 2 — Z 2 ) Y TO 2 + 02 — Z 2 



But in the sigma-model case the resolvent is not known exactly, and one must resort to an approximation such as the 
derivative expansion. Nevertheless, both cases lead to the same result (0>- „ 

At finite T, the kink case is also easy because we know the relevant part of the resolvent exactly [jlTj . Inserting (|^) 
into the contour integral expression (|f|) one finds 


(TV) 


kink 

T 


E 


2 (lM)2 sinfl 


n =0 ((2 n + l) 2 + (^) 2 ) v / (2n+l) 2 cos 2 0 + (^) 2 


(9) 


Note that (TV)|! nk is temperature dependent, but it is still topological in the sense that it only depends on the back¬ 
ground field through its asymptotic value 9 = arctan[02/m]. Also, even though the expression (^) looks complicated, 
it smoothly reduces to ([71) as T —> 0. It is also interesting that the angular nature of 9 is more obvious at finite T 
than at zero T. 
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At finite T, the sigma-model case is much more complicated, because there is no simple closed formula for the 
relevant part of the resolvent. In ||lo|| the derivative expansion was used, and it was shown that for T <C m the 
derivative expansion can be resummed to all orders, leading to the following simple expression 


(N)p m& = (V) s 0 igma 



dxe~ mt3 sinh 



+ ... 


( 10 ) 


The correction term is temperature dependent and nontopological, as it is sensitive not just to the asymptotic values 
of 9(x), but also to the details of the shape of 9{x). This nontopological term has a simple physical interpretation as 
follows [jn|. For the sigma model case we can write the interaction term in (||) as mip (cos 9 + *75 sin 9)if> = mip e llh ®ip. 
Then make the local chiral rotation ip = e~' ll5e ^ 2 ip to obtain 


_ _ _ Qt 

C = iipijfop — rmpip + ip'j 0 — ip 


( 11 ) 


which describes a fermion of mass m in an inhomogeneous electric field E(x) = —^9"{x). In this language the zero T 
induced fermion number (Jt[) can be understood as the vacuum polarization, due to the inhomogeneous electric field 
background, of virtual vacuum e + e _ dipoles. This effect is independent of the temperature as the virtual vacuum 
dipoles do not live long enough to thermalize. However, at finite T there is another effect of the background electric 
field. The real particles and antiparticles respond to this electric field in a way that can be described by linear response 
theory. In the derivative expansion limit, where 9’ <C m, the leading effect of the electric field background can be 
represented as a local chemical potential p(x) = —^9'(x). With local Fermi distributions 


f±( x j k) 


1 

e /3(VP+ij?T(i(i)) _|_ \ 


( 12 ) 


the linear response charge density is p(x) = f H f(x, k), where / = /+ - /_. For T < m, 


P(x) ~ -\[^e- m/T (J^ (13) 

in perfect agreement with the non-topological term in 

So, while at zero T, the kink and sigma-model cases each have induced fermion number given by the same expression 
(j^), at nonzero T they have very different induced fermion number, given by (J)[) and ( |l0[ ) respectively. In each case 
(N)t is temperature dependent, but in the kink case it is topological, while in the sigma-model case it is nontopological. 
The physical reason for this difference is explained in the next section. 


IV. TOPOLOGICAL VERSUS NONTOPOLOGICAL 

The straightforward mathematical reason why (N)t is topological for the kink case is that the odd part of the 
spectral function is itself topological, as shown by the Callias result (||). But a much more physical understanding 
can be obtained by considering the symmetries of the fermionic spectrum [ pd| . And this perspective generalizes 
immediately to higher dimensional models. 

Begin with the simple trigonometric identity, tanh(4p) = 1 — 2 n(E), where n{E) = l/(e^ E + 1) is the Fermi-Dirac 
distribution function. Then the finite T fermion number (^|) separates as 

/ OO 

dEa(E)sgn{E)n{\E\) (14) 

-OO 

The first term is just the zero temperature contribution, which can be expressed in terms of the spectral asymmetry, 
and is known to be topological. But the second term is highly sensitive to the details of the spectrum, because of 
the presence of the Fermi weighting factor n(\E\). Thus, the only way this term can be topological is if the odd 
part of cr(E ) is itself topological. In terms of the spectrum, this can only happen if the spectrum is symmetric under 
E —> — E . Remarkably, this does sometimes happen, but only for extremely special backgrounds. One example of 
such a special background is the one-dim kink background discussed in the previous section. Indeed, for any kink 
background (with 9 positive), the spectrum of the Dirac Hamiltonian always has the form shown in Fig. 1. 
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Paired states 



FIG. 1. The general structure of the Dirac spectrum for a kink background. 


There are continuum states for \E\ > 


m 2 + cf) 2 = m sec 9. 


In addition, there is always a bound state at E = 


—sign (6)m. There may or may not be additional bound states with m < \E\ < msec9. But if these additional 
bound states are present, they necessarily occur in ±E pairs, because of the quantum mechanical SUSY of the Dirac 
Hamiltonian in the kink case Q. Thus, the contributions of these paired bound states to the integral in (|l4|) cancel 


in pairs, while the unpaired bound state at E = — sign(0) m leads to a contribution —sign (9)n{m). The quantum 
mechanical SUSY of the Dirac Hamiltonian is likewise the key to the exact Callias index theorem result jui| for the 
odd part of the spectral function [which is equivalent to knowing the even part of the resolvent in (||)]. This means 
that the remaining integral over the continuum states can be expressed as an integral of the odd part of the spectral 
function, weighted by the Fermi factor n(E), over the positive continuum beginning at the threshold energy msec9. 
Putting all this together, we obtain 


(AOt = (iV) kmk - n(m) + 


dE 


2m 2 ta,n9n(E) 


18 7r (E 2 — m 2 )\/ E 2 — m 2 sec 2 9 


(15) 


One can check that this is, in fact, the Sommerfeld-Watson transform of the summation expression found before in 



-m m 



FIG. 2. The structure of the Dirac spectrum for a sigma model background. 


On the other hand, for the sigma-model case the fermion spectrum has no obvious symmetry. In this case, the 
Dirac spectrum has the form shown in Fig. 2. The only general thing we can say is that the continuum thresholds 
are at E = ±m. There may or may not be one or more bound states for \E\ < m. The existence of (and the precise 
location of) bound states is highly sensitive to the actual shape of the 9{x). As there is no special symmetry in the 
spectrum, there is nothing general that can be said about the odd part of the spectral function. Nevertheless, the first 
term in (|14|), which is (N) o, yields a topological result since the integral can be written, using Levinson’s theorem, in 
terms of the phase shift at infinity, which is topological ||[7|,[ij . This application of Levinson’s theorem doesn’t work 
for the T dependent second term in (0 because of the Fermi factor. This is another way to understand why the 
finite temperature induced fermion number is generically nontopological. 

These arguments based on the symmetry of the fermion spectrum (in the given background) clearly generalize to 
any dimension of spacetime. So, we conclude that the generic situation is that the finite T induced fermion number 
is non-topological; but, if the background is such that the fermion spectrum is highly symmetric, then the induced 
fermion number may be topological even though it is temperature dependent. 


V. THERMAL FLUCTUATIONS OF FERMION NUMBER 

Fluctuations in the induced fermion number can be analyzed in an analogous manner [ |i~T| : 

(A1V) 2 = ( N 2 > - ( N ) 2 = 11 °° dEa(E) sech 2 (^j 
The trigonometric identity, |sech 2 (4^) = n(E)( 1 — n(E)) 1 re-expresses this as 


(16) 
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cLEcr(E) n(E) (1 — n{E )) 


(17) 


/ OO 

-CO 

which is simply a sum over the fermion spectrum of n(l — n), as it should be for the thermal occupation of independent 
single-particle states Q. From it is clear that the fluctuations are always non-topologicah as they are sensitive to 
the details of the fermion spectrum. Another way to see this is to note that the expressions ( [I fi[|l 7| i require knowledge 
of the even part of the spectral function, which is not accessible to an index theorem. So the fluctuations must 
typically be calculated approximately; however, for a special class of kink backgrounds, </> 2(2) = 4> tanh(</>a:/j), where 
j is an integer, it is possible to compute {AN) 1 2 3 4 5 6 7 8 9 10 exactly, and one sees the explicit (nontopological) dependence on the 
scale factor j |ll|| . 

These expressions also make it clear why the fluctuations vanish at zero T, because the sech 2 {(3E/2) factors vanish 
exponentially fast as T —> 0. Thus, the induced fermion number is a sharp observable at T = 0, in agreement with 
previous results [p~8[ 7 but at nonzero T it is not a sharp observable, simply because the thermal expectation value 
mixes states other just the ground state. Another example of non-sharp fractional fermion number (but not in the 
context of temperature dependence) has been discussed recently in the context of liquid helium bubbles p9[ . 


VI. CONCLUSIONS 


To conclude, the finite temperature induced fermion number (N)t is generically nontopological, and is not a sharp 
observable. The induced fermion number naturally splits into a topological temperature-independent piece that 
represents the effect of vacuum polarization on the Dirac sea, and a nontopological temperature-dependent piece that 
represents the thermal population of the available states in the fermion spectrum, weighted with the appropriate 
Fermi-Dirac factors. As the temperature approaches zero, the nontopological terms vanish exponentially fast, and we 
regain smoothly the familiar topological results at zero temperature. But at nonzero temperature, the nontopological 
contribution is more sensitive to the details of the spectrum, and so is generically nontopological. This argument 
applies in any dimension. 

Consider, for example, for a chiral SU{2) sigma model with a Skyrme background in 3+1 dim., 


in’ll) (ttq + * 757 ? • t) if) 

'1 


mi/> 


2 (f? + 5 f ) + 2 (3 





where 7Tq + if 2 = 1 and 7To + *7? • r = g , which is the analogue of the 1+1 dim. sigma-nrodel case analyzed in this paper. 
From numerical work pl[ |, the Dirac energy spectrum of the fermions is not symmetric, and the energy of a possible 
bound state (and indeed the number of bound states) is highly sensitive to the details of the shape of the radial 
hedgehog field. This shows that in this case the finite temperature induced fermion number is nontopological. Thus, 
for a (large) Skyrme background in 3+1 dim., the finite temperature fermion number is not simply the (topological) 
winding number of the Skyrme field, but a much more complicated nontopological object. 
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